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Abstract 

Physically relevant gauge and gravitational theories can be seen as special 
members of hierarchies of more elaborate systems. The Yang-Mills (YM) sys- 
tem is the first member of a hierarchy of Lagrangians which we will index by 
pi , and the Einstein-Hilbert (EH) system of general relativity is the first mem- 
ber of another hierarchy which we index by p2- In this paper, we study the 
classical equations of the pi = 1,2 YM hierarchy considered in the background 
of special geometries (Schwarzschild, deSitter, anti-deSitter) of the p2 = 1,2,3 
EH hierarchy. Solutions are obtained in various dimensions and lead to several 
examples of non-self-dual YM fields. When pi = p2 self-dual solutions exist 
in addition. Their action is equal to the Chern-Pontryagin charge and can be 
compared with that of the non-self-dual solutions. 
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I Introduction 



The main aim of this work is the study of non-selfdual Yang-Mills (YM) fields in 
d > 4 dimensions on fixed gravitational backgrounds in (i-dimensions, extending the 
work of 0, ^ in 4-dimensions. As in |1|, we restrict to Schwarzschild and 

deSitter metrics. By YM fields here we mean the solutions to the hierarchy of p- 
YM systems |^ , whose p = 1 member is the usual YM system, and the generic p-th 
member involves the curvature 2p-form in lieu of the usual YM 2-form curvature. 

Here we present nonselfdual solutions in curved backgrounds. Such solutions in 
flat spaces would be particularly interesting. Even the complex one in dimension 4, 
which was obtained p|, |^ exploiting the conformal flatness of the deSitter metric in 
this case, is worthwhile. (See the remarks and the references in concerning the 
relevance of complex saddle points.) 

We have a definite reason for considering the hierarchy of YM systems rather 
than restricting to the usual {p = 1) member only, other than the fact that it is quite 
natural to do so in dimensions d > 4. As highlighted in the work of Refs. ^ |^, 
in c? = 4 dimensions and with Euclidean signature, the selfdual {p = 1) YM field can 
be constructed from the double-selfdual gravitational field by constructing the SU{2) 
YM connection from the corresponding gravitational spin-connection [^. In this case, 
it is known that the gravitational field equations are automatically satisfied, and so are 
the YM equations by virtue of the selfduality0. Clearly then, the YM action will be 
equal to the Chern-Pontryagin (C-P) charge. Thus, if one finds non-selfdual solutions 
to the YM equations in the said double-seldual gravitational background, the value 
of the action will differ from that of the self-dual YM fields, and it is interesting 
to compare it with the value of the C-P charge originating from the double-selfdual 
metric of the previous case [||, |]. For the usual (p = 1) case a selfdual YM field can 
thus be related to a double-selfdual gravitational metric only for d = 4, and it is our 
desire to carry it out in > 4, as a subsidiary motive of our investigations, that leads 
us to consider the hierarchy of YM systems. 

The main interest of the present work remains the investigation of non-selfdual 
YM fields. In fiat 4 dimensional Euclidean space no explicit real nonselfdual solutions 
are known, cf. the explicit complex non-selfdual solutions in 0] for gauge group 
SU{2). (Implementing instantons and antiinstantons in commuting subgroups of 
SU{N) for sufficiently large one can evidently obtain nonselfdual solutions, but 
here we consider specifically SU{2).) For SU{2) only selfdual solutions are known, 
the value of the YM action pertaining to these being equal to the Chern-Pontryagin 
(C-P) charge. It was found in |l| however that on fixed Schwarzschild and deSitter 
backgrounds, there were nonselfdual solutions, the values of whose actions differed 
from that of the C-P charge. In particular, in the case of Schwarzschild background, 
the value of the nonselfdual action, which turned out to be real, was slightly smaller 
than that of the selfdual action. Thus, the study of YM fields on fixed backgrounds 
is of general interest in the context of nonselfdual fields, and in the special case(s) 
where selfdual solutions on the same background also exist, then it is of particular 
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interest to compare the value of the nonselfdual action with that of the selfdual one. 
The reason for considering the hierarchy of YM models is precisely because the p- 
th member of these, on Ap dimensions, does support selfdual solutions. Indeed, the 
situation in 4 dimensions, where the YM field constructed from the double-selfdual 
spin-connectionon is automatically (single) selfdual, occurs in all 4p dimensions for 
the p-th member of the YM hierarchy on a background whose 2p-form Riemann tensor 
is double-selfdual. We shall refer these 4p-dimensional gravitational systems as the 
hierarchy of Einstein-Hilbert (EH) systems, or, generalised EH systems |p . 

The hierarchy of (gravitational) EH systems in all even dimensions was previously 
studied in Refs. and p|, and more recently in 0. In |0] it was shown that in 4p 
dimensions, the 2p-form Riemann tensors of the deSitter and Fubini-Study metrics 
were double-selfdual, and in |^ this was done for Schwarzschild like metrics. In 
it was shown that the YM connections constructed from the corresponding spin- 
connections in all even dimensions yielded 2A^-form YM curvatures which satisfied 
the (single-) self duality conditions of the YM hierarchy This construction is pos- 
sible in all even dimensions only when the curved space is a compact coset space, but 
in general can also be carried out if we restrict the dimensions to Ap. In the present 
work, where YM fields (both selfdual and nonselfdual ) on Schwarzschild backgrounds 
are considered, we restrict to Ap dimensions, and we label the EH systems also with 
the label p. 

Since we will consider the p-th members of both YM and EH hierarchies in all 
d dimensions, with d ^ 4p for the non-selfdual cases, we will henceforth label the 
members of the YM hierarchy with pi and the members of the EH hierarchy with 
P2- For a discussion of the YM and EH hierarchies, we refer to ^ and p], but here 
we simply give the definitions for the Lagrangians of these systems for the particuar 
examples that will be employed in the present work. 

Thus, concerning the YM systems, the pi = 1 and the pi = 2 systems read 
respectively 

5(1) = TrF'^'^F,,, (1) 
5(2) = TiF^-'P'^F^.^p,. (2) 

with the following definitions for the field strenghts (the potentials are antihermitian, 
the gauge group will be specified later) 

F^, = d^A,-d,A^+[A^,A,] , (3) 

F^upcT — -^pcr] } (4) 

and the square bracket on the indices [i^pcr] implies cyclic symmetry. The final normal- 
isations of both (0) and (^ will be fixed such that the value of the action pertaining 
to the spherically symmetric self-dual solutions on Schwarzschild background (with 
Euclidean signature), is set equal to one. 

Concerning the EH systems, these are the p2 = 1, P2 = 2 and the P2 = 3 systems 
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p _ ui..ti4iyi..Ud-4,f,ni na-icmi..ii4ni..nd-2 pmim2 jDmzmA, (a\ 

f _ ii\ii2..ii(,vx.Ma.-ii n^ na-ec"^i--me,ni..nd-2 pmim2 pm3m4 pmsmg fn\ 

'^(3) — ^ '^i^l -^MlM2 -^^3/^3 -^M5M6 ' 

in d dimensions. The corresponding (generalised) Einstein equations are given in 

In the self-dual cases, with pi = P2 = d/A, the stress tensor due to the YM fields 
vanishes identically so that the (generalised) Ricci scalars (|§)-(0) vanish for these field 
configurations. Thus the action of the YM fields equals the Chern-Pontryagin charge. 
To calculate the action integral in (static) Schwarzschild backgrounds in such cases, 
one has to integrate the time variable over one period. This is the period associated 
with the desingularisation of the Schwarzschild metric by introducing Kruskal like 
coordinates in the case of Euclidean signature. For the arbitrary dimensional case, 
and for the dynamics determined by the p2 EH system, this period was calculated 
to be 

where /C is a parameter [rU| in the (hierarchy of) Schwarzschild metric(s)given in H]. 



In our calculations of the actions below we will suppress the factor (||) contributed by 
the (Euclidean) time integration, since we are only interested in the relative values of 
the selfdual and non-selfdual actions for any given EH-YM system. 

In m, very simple nonselfdual solutions for the d = 4 (usual) pi = 1 YM fields 
in Schwarzschild and deSitter backgrounds (both of the p2 = 1 member of the EH 
hierarchy) in four dimensions were presented. They were further discussed in Refs. 
and 0. Here they are generalized to all dimensions d > 4, with Lorentz or Euclidean 
signature and {d — 1) spatial dimensions, for both pi = 1 and pi = 2 members 
of the YM hierarchy, on the backgrounds of the p2 = 1,2,3 members of the EH 
hierarchy. The corresponding constructions for members of the pi > 3 on various p2 
backgrounds can be given systematically. As in the 4-dimensional pi = 1 YM case [|I| 
on p2 = 1 Schwarzschild and deSitter backgrounds, where these solutions are found 
in closed form, it turns out that in the d-dimensional {d > 4) pi = 1 YM case on 
these P2 = 1 backgrounds, non-selfdual solutions can also be constructed analytically 
in closed form using the same procedure. On p2 > 2 backgrounds however, non- 
selfdual solutions of the pi = 1 YM system could only be eveluated by numerical 
construction. In the pi = 2 (as well as all pi > 2) YM case(s) likewise, the non- 
selfdual solutions on p2 EH background could be constructed only using numerical 
integrations irrespective of p2- These are the main results of the present work and are 
presented in Section 2. The corresponding constructions for members of the p > 3 
can be given systematically. 

For the pi = 1 and pi = 2 systems in c? = 4 and d = 8 dimensions respectively 
(with Euclidean signature), the selfdual solutions are considered in detail for P2 = 1 
and p2 = 2 Schwarzschild metrics, respectively. These results are needed for the com- 
parison of the Euclidean actions of the selfdual and nonselfdual solutions, which we 
carry out for these cases. The corresponding analysis with deSitter backgrounds is 
not carried out in detail since in that case there are no (real) non-selfdual solutions. 
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Again, extension to the generic (pi = P2)-th self dual cases can be carried out sys- 
tematically. These results are presented in Section 3. In Section 4, we give a short 
summary and discussion. 

Solutions on the anti-deSitter backgrounds of p2 = p EH, to the pi = p YM in 
d = 4p dimensions, are studied in the Appendix. 

II Non— selfdual solutions on fixed backgrounds 

This Section is subdivided in three Subsections. In the first one, we state our Ansatz 
for the YM fields on fixed Schwarzschild and deSitter curved backgrounds using the 
Kerr-Schild parametrisation of these metrics and give the Euler-Lagrange equations 
whose non-selfdual solutions we seek. In the second and third Subsections, we present 
the non-selfdual solutions on Schwarzschild and deSitter backgrounds, respectively. 

Since all the actual calculations involved in this, and the next Section, involve 
only Euler-Lagrange equations and selfduality conditions of members of the YM 
hirerchy, we will not be engaged in a description of the hierarchy of Einstein-Hilbert 
(EH) systems. The only information we need here in this respect are the actual 
functions parametrising the Schwarzschild and deSitter metrics pertaining to these 
hierarchies, which we have stated when needed below in (|16D and i\i7\j . For details of 
their derivations, we refer to Similarly for the hierarchy of YM systems we simply 
state their field equations (|32D-(p3D and, refer to for their general details and to 
§ for their properties relative to the EH hierarchy. 

II. 1 Ansatze and YM equations 

There are several ingredients needed for making the spherically symmetric Ansatz. 
One is the definition of spin matrices, or the gamma matrix representations of the 
generators of SO{d) (the choice of this gauge group is dictated by our requirement of 
spherical symmetry in d and d — 1 dimensions) . Another ingredient is the parametri- 
sation of the components of the metric for the (fixed) curved space on which the YM 
fields are defined. The ansatz for the gauge potentials has then to be specified and 
finally the equations derived. We deal with these items in that order below. 

II. 1.1 Gauge group and representation 

Even though our dominant motivation here is the study of non-selfdual solutions, 
let us start with the case of selfdual solutions which we shall consider in the next 
Section for comparison with the results of this Section. This puts a restriction on 
the 2^^/^ X 2'^/^-dimensional representations of the SO{d) matrices in spacetime with 
d-dimension. The representations of the gauge groups of the YM fields pertaining 
to the p-th member of the YM hierarchy are chosen such that in Ap dimensions, 
there exist selfdual solutions on R''^ namely that the gauge group is represented 
by 2('^-2)/2 X 2('='-2)/2 = 22p-i X 22p-i (left or right) chiral S0{4:p) matrices, denoted 
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by S0±{4:pi). This choice of gauge group representation also makes it possible to 
construct single-self dual pi = p -YM fields oiap2 = p -EH backgrounds @| . Examples 
of the latter are known [0, ||, |] , and of these the deSitter and Schwarzschild 
will concern us in this work. 

The representations of the algebra of S0±{4p) employed in our Ansatze are de- 
noted by 

Safe = -^S[,Sb] = -i(Il±r4p+i)[r„,r,], (a,6= l,2,..,rf) (9) 

where the square brackets on the 4p component indices [ab] imply antisymmetrisation, 
Ta denote the d = 4p -dimensional Gamma matrices and r4p+i is the corresponding 
chiral matrix. 

As we shall be concerned with p = 2-YM systems in detail in this and the next 
Section, it is convenient to define the totally antisymmetric 2p-form tensor-spinor 
matrix 

^abcd = {Sa[fe,Scd]} (10) 

using the same notation as in (^). In addition, we state an identity which will be 
useful in the next Section, 

{^ab, ^cd} = -^abcd — -{^ac^bd — SadSbc)'^ ■ (H) 

Let US now relax the constraints necessitated by the requirement of YM systems to 
support self dual fields, thus allowing d dimensional SO{d) systems, without requiring 
that d = Ap. For even d this could allow the assignment of chirally symmetric 2^^/^ x 
2d/2 representations for the spin matrices. These are not the representations we assign 
to the even d dimensional gauge algebras, but restrict to the chirally asymmetric 
(left /right) 2'^'^"^^/^ x 2'^'^"^^/^ spin matrices defined by (^. The reason is that these 
are the relevant representations for the selfdual cases, which we wish to compare 
eventually against the non-selfdual cases under consideration. For the static field 
configurations (deSitter or Schwarzschild) we will be concerned with, the spherical 
symmetry will be imposed in the odd, d — 1, dimensions, when there there exits no 
chiral matrix. The representations of the spin matrices in this case are defined by the 
2(d-2)/2 ^ 2*^'^"^)/^-dimensional SO{d— 1) matrices Sj^, with = 1, 2, ..,d~ 1) given 
by (I), or equivalently by 

s^, = -J[r.,r,] (12) 

in terms of the (d — 2) Gamma matrices Fi, .., rd_2 (with dimension 2*^°'^^^/^ x 
2(<i-2)/2-j^ supplemented by their chiral matrix F^-i. In what follows, the precise 
dimensionality of the representations of Sj^ will not matter except in the (canonical) 
dimensions where selfdual YM fields are supported, inosfar as the value of the action 
densities depend on these. Otherwise their only important feature will be the fact 
that they satisfy the algebra of SO{d) or SO{d — 1) respectively. 
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II.1.2 The metric 

Next, we give tlie Kerr-Scliild parametrisation of the background metric, which will 
be employed in the Ansatz. We consider spherically symmetric, static metrics for 
dimensions d > 4 

9^.. = V^u + l^.lu, g''" = 7]^^" -Fr (13) 

where 

^700 = = -1 , Vij = V'^ = ^ij 

and 

n-'V^, = FFg,, = . (14) 
For static spherical symmetry Iq is a function of r only, where 

d~l 
2 \ ^ 2 

i=l 

and 

X ' 

li = lo— = loXi, {i = l,2,...,d-l) 
r 

satisfying 

-ll + T.l" = (15) 

i 

The Schwarzschild metric in d dimensions pertaining to the usual {p2 = 1) EH 
system, was given by Myers and Perry |Tl|. The correspoding result in d dimensions 
pertaining to the generic p EH system was given recently in Ref . |^ . For the generic 
case we have 

ll = Cr-(^^) , (C > 0) (16) 
and for deSitter metric, for all d pertaining to all p EH systems, 

ll = Ar^ (A > 0) . (17) 
The standard form in spherical coordinates is given by 

ds^ = -Ndt^ + N-^dr"^ + r'^dn^d-2) , (18) 
where dQ(^d-2) is the line element on the unit {d — 2)-sphere and 

N = {1- ll) (19) 
The coordinate transformation relating (0) and (|T8|) , namely 

/dr 
j^-r (20) 

does not affect our particularly simple ansatz for the gauge potentials to follow (with 

At = Ar = 0). 

After constructing the solutions using ([T^) the passage to Euclidean signature is 
best considered ( rather than introducing imaginary /q) by directly starting from (|T^, 
leading to 

ds"^ = Ndt^ + N-^dr^ + r'^dVL^d-2) (21) 
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II. 1.3 The gauge potentials 

The ansatz for the gauge potentials is 



Ao = (22) 
A, = (ir(r) - l) (^ = l,2,...,rf-l) (23) 



with Ejj defined according to (||), or equivalently (|12]), to be the spinor representations 
of S0{d—1). The corresponding YM curvature (|]) has then the following components 

Fo, = (24) 
Fij = Vi T^ij + V2 X[i Xk (25) 

where as before, the square brackets on the subscripts [ij] imply antisymmetry, and, 

V, = -r-\K^-l) (26) 
V2 = -r-^[rK' - {K'^ ~ 1)] . (27) 



Now using (|T3|) we have 



with 



= llK'T.,kXk (28) 
= g^^g^'^F^, 

= Wii:,^ + W2X[iEj]kXk, (29) 



Wi = -r'\K^-l) (30) 
W2 = -r'^[rNK' - {K^ - 1)] . (31) 



II. 1.4 The equations 

The Euler-Lagrange equations of motion are expressed very simply since | (7 |= 1 for 
the KS metric. For the p = 1 and p = 2 YM systems they are, respectively 

D^F^"' = (32) 
{Fp., D^F^'-P''} = 0. (33) 

with the definitions (D,(D. 

After some straightforward but laborious calculations, substituting (p2|), (p3D , and 
(|28|), (|29|), into ( |32D and (|33|), we find first that the Gauss law equations 

= and {Fp^ , = 
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are identically satisfied, and the remaining components of ( P^j) and ( P^D yield simply 



{[Nr'^-'^K']' -{d- 3y-^K{K^ - 1)) x S^-fc (34) 



Vi {[Nr'^-''{K^ - l)K']' -{d- by-'^^KiK^ 
X 3{d-3){d-A) J^jk Xk . 



X 

(35) 



The Euler-Lagrange equations of the one dimensional subsystems of the d dimen- 
sional p = 1 and p = 2 YM systems in the static spherically symmetric background 
specified by (|TBp can readily be read off ( P^D and (|55|), as single equations for the 
functions K{r). Indeed, it is easy to read off the equations for the arbitrary p YM 
systems: 

[Nr'^-^^K^ - ly-^K']' =[d- {2p + l)]r'^-2(2p+i)ir(ir' 
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(36) 



We discussed in the Introduction, the Einstein-Hilbert (EH) hierarchy in parallel 
with the YM hierarchy. The former determines the metric on whose background 
the YM fields are studied. In the context of (|36D , it is the function N{r) that is 
determined by the dynamics of the relevant EH system. Since in this Section we 
are interested in solutions on fixed backgrounds, there is no reason to privilege any 
particular member of the EH systems for characterising the function given by ([16D , 
([T7| ) and (p^l). Using the notation introduced in Section 1, we label this function 
by the index p2 pertaining to the EH hierarchy, namely as N(^p^^, and simultaneously 
rename the index p in ( |36| ) pi. Equations (^) now are expressed more specifically as 

[Nip^y-^^'iK^ - ly-^K']' =[d- (2j9i + l)y-^^^P'+^^K{K^ - . (37) 

For use in comparing the actions of non-selfdual and selfdual p = 1 and p = 2 
YM systems in c? = 4 and d = 8 respectively in the next Section, we present the 
action densities of these two systems for the fields (p4|) , ( p5|) , (|28|) and (^), for which 
we will give the non-selfdual solutions in the following two Subsections. Because of 
the vanishing of F^i, (p^), these action densities coincide with the (Euclidean) energy 
densities. They are, for the p = 1 YM system. 
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Tr F^^F^, = Tr F^^F,, 
2'^{d- 2)[{d - 3)WiVi + 2{Wi 



W2)iVi-V2)] 



d-6 

2— 



id -2) 



2r'NK" + {d-3){K^ - I] 



and for the p = 2 YM system. 



5, 



(2) 



Tr F^^'"^F, 



Tr F^'^^Fh 



d-lO 

2— 



{d-2){d-3){d-4:) X 
xWiVi[{d - 5)WiVi + 4(W^i - 
2^(rf- 2)(rf-3)(d-4) X 



(38) 



(39) 



(40) 



W2)iV,-V2) 



X 



(K' - 1] 



Ar'NK" + {d-5){K' -1) 



(41) 
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In (^) and (PI), the numerical factors resulting from the traces of the S matri- 
ces are accounted for. We note here that varying the densities r'^~^Tr F^^Fij and 
j-'^-'^Tt F''^^^Fiijk, given by ( |39D and ( PD respectively, with respect to K{r), we obtain 
the equations { ^7\j with pi = 1 and pi = 2 respectively. This is not surprising. 

Before proceeding to integrate (|36D for p = 1 and p = 2, we make some general 
remarks. Firstly, equations (|36|) do not satisfy the Painleve criterion |12|, nonetheless. 



we find some special solutions. Perhaps the most remarkable feature of equations (pq) 
is the fact that on flat background, with = 1, there are no solutions so that the 
only non-selfdual solutions are on curved backgrounds. It is interesting to make some 
general observations at this point. We note that in the generic case, in dimensions 
c? = 2p + 1, the right hand side of (^) vanishes. In particular, for the p = 1 this 
coincides with d = 3 corresponding to the Abelian case with only E12. In what follows, 
we will restrict our attention to d > ip for each case p. Special features arising for 
d = ip + 1 will be discussed also. 

II. 2 p2 Schwarzschild backgrounds 

We now present the solutions of Eq.(|37|) for the Schwarzschild case i.e. for 



d-2p2-l _ d-2p2-l 

iV(,,) = !-(-) , C^C^^ (42) 

(where we redefined the parameter C of (0)). We will limit ourselves to those values 
of P2 for which the metric (pl|) is asymptotically fiat, i.e. to p2 < {d — l)/2. Eq.(|3 
3 to be solvec 



has to be solved on the interval [C, 00] with the boundary conditions 

-CK'(C) -id- 2p2 - l)KiC){K^iC) - 1) = , Kioo) = 1 (43) 

which arise by demanding the regularity of the solution at r = C* and the finitenes of 
the action. 

II. 2.1 Pi = I , P2 = I case 

Solutions to (|37| ) with one-node can be constructed by using the ansatz 

This leads to simple algebraic constraints on the parameters a, b 

3a + b{2d-7) + (d-l) = (45) 
a{a + b)-{d-5)b = . (46) 

which involves essentially solving only a quadratic equation. For d = 4 the old results 
in are reproduced. Of the two real solutions only the one with 

a = -2.366, b = 4.098 (47) 
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gives finite energy (action) for Lorentz (Eu- 



(the exact values follow from (53)- 
clidean) signature. 

For d = 5 one has a very special CclSG, clS IS evident from (^B]). The solution 



0, 



■(4/3) 



(4J 



leads to a divergent action since the domain of (r/C) is [l,oo] and this includes a 
zero of the denominator of K. But now one can also consider the flat limit as follows. 
Setting, for (rf = 5), 

a = d/C^, b = b/C^ (49) 



and (^61) reduce to 



3{d + b)+AC^ 
d{d + b) 





0. 



Hence as C 



one obtains 



Substituting 



0, setting (5 being an arbitrary real number) 

b = 5^ 



m 



— a 

{K - 1) : 

(g§)and 

ei234 = 1 , 



(50) 
(51) 



(52) 
(53) 



r2 + 52 

) it is seen that for the convention, say, 
-S34 (cychc) 



for the chirally projected 2x2 S0{4) generators one obtains the famous BPST 
selfdual solution in d = 4 as the static limit in c? = 5 via our limiting process. 
Another convention gives the antiselfdual form. 

/^From d = 6 onwards the solutions become complex. The corresponding finite 
complex action ,or energy, will be obtained in the following section. The exact values 
can be obtained, for any d, immediately from (^3p-(^B|). Some numerical values, 
giving a direct idea of variation with d, are given below. Both upper or both lower 
signs are to be taken. For 



d = 


6 : 


a = 


0.500 ±zl. 500, 


b = -1.300 TzO.900 


d = 


7 : 


a = 


0±il.732, 


b = 0.857^^0.742 


d = 


8 : 


a = 


-0.167 ±a.863. 


b = -0.722 ^i0.621 


d = 


9 : 


a = 


-0.250 ±a.984. 


b = -0.659 TiO.541 


d = 


10 : 


a = 


= -0.300 ±i2.100. 


b= -0.623 TiO.485. 



We now come back to the case d = 4 for which we were able to further construct 
numerically a solution with two-node. It has in particular 



K{1) ^ 0.045 , K{2.35) ^ 0.0 , K{54:. 



0.0 



(54) 
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and the two profiles of K are plotted on Fig. 1. In this figure, it is seen that the two 
solutions for K{r) cross the value iiT = at nearly the same value r 2.36 (although 
numerically different) and that the two-node solution reaches its asymptotic value 
K = 1 very slowly. 

The occurence of such a couple of solutions suggests that the action functional 
(see Eq. (|55|) below) admits an infinite series of extrema indexed by the number of 



nodes of the function K{r). This is reminiscent to the series of Bartnik-McKinnon [|T3 
solutions in Einstein- Yang-Mills theory. Similar series of particle-like solutions have 
also been discovered to the Einstein- Yang-Mills-Higgs equations [Q. For a review, 
see 
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The construction of the (eventual) additional solutions in our case would demand 
more numerical analysis and lies outside the scope of the present work. 

The different solutions can be characterized by their action. Up to trivial factors 
this is determined by the following radial integral 

/oo 
dxx''-^[2r^Ni^,){K'f + {d-?,){K^ -If] . (55) 

For the two solutions of the case ci = 4, we find 

jf ^ ^ 0.959 , jf ^ ^ 0.992 (56) 

respectively or the one-node and two-node solutions. 

No solutions other could be constructed numerically for the equation ( pTj ) with 
p = 1 and d > A. This is probably because, if any, they are not real. 



II.2.2 pi = 2, p2 = 2, 3, 4 case 

In the case pi = 2, we could not find a self consistent ansatz like (^4|); however, we 
manage to construct a few real solutions numerically. In analogy with the case pi = 1 
above, where the only real solution was in d = Api = 4 dimensions, here we expect to 
find real solutions in d = 4pi = 8 dimensions. This indeed turns out to be the case. 
In addition, it is possible to employ the three members labeled by p2 = 1? 2, 3 of the 
EH hierarchy to specify the fixed background. 

The non trivial factor of the action integral reads in this case 

/oo 
dx x^-i°[4x2iV(p,)(ir')'(^' - 1)' + id- 5)iK' - If] , (57) 

We were able to construct numerically a solution where K{r) develops one node 
and another one where K{r) develops two nodes. For the three cases p2 = 1, 2, 3, the 
numerical evaluation of the action/energy integrals (|57|) leads to 

J92 = 1 1^ ^ 2.61 ^ 2.95 (58) 

P2 = 2 ^ 2.38 1^ ^ 2.90 (59) 

P2 = 3 1^ ^ 2.17 iP ^ 2.81 (60) 
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Again, the existence of two-node solutions suggests that the functional (0) admits 
an infinite series of extrema indexed by the number of nodes of the function K{r). 

Inspection of (pSj), ( |59| ) and ( |60D reveals that the action/energy integral /(s) of the 
Pi-th member of the YM hierarchy on the fixed background of the p2-th member of 
the EH hierarchy, decreases with increasing p2. This is true for both one-node and 
two-node solutions, and we expect it is a general feature for such sytems. 

The profiles of the solutions are presented in Fig. 2. Again in the case p2 = 1 we 
remark that the two solutions cross K = for nearly the same value (r ^ 2.13) of 
the radial variable. This phenomenon is less and less true for increasing p2. 

Our numerical analysis of the equations for d > 8 has lead to no other solution. 

Superpositions (by means of linear combinations) of the lagrangians with different 
values of pi and/or of p2 could be considered as well, leading to a many-parameter 
differential equation but we have not considered such possibilities. 



II. 3 deSitter backgrounds 

In this case, there is one single background function N{r) for all members p2 of the 
EH hierarchy, in all dimensions d. The relevant equation to solve (^7D with 



N={l-Ar^) (61) 

(see ( p!?!) and (p!9D). We therefore consider only the different members of the YM 
hierarchy and label them with pi = p throughout this section. 

The equation has to be solved on the interval r G [0, A^^/^] with the boundary 
conditions 

K{0) = 1 , 2AK'{A) + {d-2p-l)K{A){K\A)-l)=0. (62) 

where for brevity we have used A = A^^/^. 

Again, in practice there is a marked difference between the p = 1 case and all 
others. It turns out that for p = I, very simple solutions can be found in closed form, 
while for all other cases the solutions can be constructed only numerically. 



II. 3.1 p = 1 case 

In this case, solution can be constructed algebraically by using the ansatz 

1 + aAr^ , , 

Substitution of this form into the equation leads to the following conditions for a, b : 

{d - 3)a{a + b) + 2{d ~ 5)b = (64) 
3{d - 3)a ~ {d - ll)b + 2{d - 1) = 0. (65) 
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Exact solutions can then be obtained by solving a quadratic conditions in a or b. 
Approximate numerical values are presented below. For d = 4 the old results |l|, 
are reproduced with 

a = ±il.732, 6= -0.857 TiO. 742. (66) 

Unlike in the Schwarzschild case, the d = 4 dimensional solution of (|36D is complex. 
For the special case d = 5 the eqations reduce to 

a{a + b) =0 (67) 

3(a + 6)+4 = (68) 

whose only consistent solution, 

a = 0,6 = -(4/3) , (69) 

leads to divergence in K since one considers the domain < Ar^ < 1. 

/^From d = 6 onwards (exhibiting a behaviour complemetary to the Schwarzschild 
case) the solutions become real. The solutions read, with e = ±1, 



-{d - 3)(rf - 4) + e^{d - 3)(5d-23) ^^^^ 



(d-3)(rf-5) 



-{d^ -9d + 26) + 3e^{d-3){5d-23) ^^^^ 



{d-ll){d-5) 

Only the solution corresponding to e = —1 leads to a regular solution on r G [0, A] 
and for 5 < d < 11. The single zero of the function K{r) reads immediately from the 
ansatz and the fact that a_ is negative. 

When solving numerically the equations for floating values of d for d 5 and 
d ^ 11, we got evidences that the solution is running into problems. 

The parameters a, b of the regular solutions together with the value of the ac- 
tion (in fact of integral /(rf) (55) now taken on r G [0,A~2]) for the different d are 
summarized in Table 1. 



Table 1: The values a, b and the action for the solutions in the case p=l. 



d 


a 


b 


I 


6 


-3.527 


4.349 


2.506 


7 


-2.366 


4.098 


1.917 


8 


-1.948 


5.073 


1.649 


9 


-1.729 


7.558 


1.496 


10 


-1.593 


15.448 


1.400 



Again unlike in the Schwarzschild case, we did not find any solutions with two or 
more nodes, despite our (numerical) efforts to do so. While this may signal the fact 
such solutions on deSitter backgrounds do not exist, this is not necessarily case. In 
the latter case, it would be a challenge to find the multinode solutions. 
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II.3.2 p = 2 case 



In analogy with the p = 1 case above, we would expect that there exist no real 
solutions to (|36D in d = 4p = 8. Similarly, we would expect that the solution in 
d = 4p + 1 = 9 would have divergent action/energy integral /(g). Accordingly, we 
would expect to find (real) solutions for > 4p + 2 = 10. This turns out to be the 
case. Surprisingly, the solution for ci = 11 turns out to be of the closed form (pBD, 
with 

^(^) - ■ 

This is in contrast to the p = 2 Schwarzschild cases where all solutions were con- 
structed numerically. All the other solutions can be constructed only numerically. 

We found numerical solutions with one node of the function K{r) for 10 < c? < 16, 
recovering the solution (p2| ) numerically. The numerical approximation of their energy 
and of the position tq of the node of K{r) are given in Table 2. We do not exhibit 
the profiles of the functions K in these cases because they are very close to the profile 
for the (analytically known) pi = 1 solution. 



Table 2: The values of the node of the solutions and the action for the solutions in 
the case p=2 



d 


VAro 


I 


10 


0.61 


3.56 


11 


1/V2 


2.67 


12 


0.76 


2.23 


13 


0.79 


1.96 


14 


0.82 


1.79 


15 


0.84 


1.66 


16 


0.85 


1.57 



To summarize, we found one-node solutions on deSitter background in 4p + 2 < 
(i < 8p + 2 for the p = 1 YM systems, and 4p + 2 < d < 8p for the p = 2 YM systems. 



Ill Selfdual YM on double-selfdual EH 

In Ref. P], the result of Ref. P] for the p = 1 case has been extended to the arbi- 
trary p case. This states that in 4p dimensions, the {pi =)p-th member of the YM 
hirearchy on the double-selfdual background of the {p2 =)p-th member of the EH 
hierarchy is selfdual. Unlike the non-selfdual solutions studied in the previous Section 
which satisfied the YM equations (|37D on a fixed background, these solutions satisfy 
the full gravitational-gauge field system taking into account the backreaction of the 
gravitational system on the YM field |1T6|. 
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Our aim in this Section is to calculate the action densities of the p = 1 and p = 2 
YM systems in d = 4 and d = 8 respectively, for the purpose of comparing their 
actions with those of non-selfdual solutions of these systems found in the previous 
Section. 

The double-selfdual 2p-form Riemann curvature in d = 4p, whose metric auto- 
matically satisfies the variational equations of the p-th EH system yields 
the 2p-form (single) selfdual YM curvature in the chiral representation so±{4p) of 
S0{4p) that satisfies the variational equations of the p-th YM system. Following 
Refs. §, H this YM curvature is given by the Riemann tensor as follows 

Ff,u = -- R^t ^ab (73) 

where the Greek letters fi, v are the coordinate indices and the early Latin letters a, h 
the frame indices, both running over (1, 2, 3, 4p). Renaming 4p as 0, the coordinate 
index [i runs over (i, 0) and the frame index a over (m, 0), where the Latin (coordinate) 
indices and (frame) indices m, n run over l,2,..,4p— 1. in ( [73D is defined 

by (1). 

While both Schwarzschild and deSitter metrics result in double-selfdual 2p-form 
Riemann tensors, it is only in the background of Schwarzschild that the YM fields have 
real non-selfdual solutions, which were found in the perevious Section. Since selfdual 
YM fields on double-selfdual backgrounds are real, it follows that for the purpose 
of comaring the latter with the former, only the Schwarzschild case is relevant to 
the work of this Section. The EH systems supporting double-selfdual solutions in 4p 
dimensions has been studied in Ref. [0, which we exploit here. 

Using the the Kerr-Schild parametrisation of the last Section, but now following 
the convention of Ref. |Q| with C replaced by 2C and L = Iq, the components of the 
Riemann tensor can be readily calculated 

RT = '-^^Pl^-jiiL'-'-^)+CLL']£,6p-^ (74) 

CT' [' CT T' 

^ il-CL) 5T + C[{L" --) + ]x,x"' (75) 

RT = ^— ^[^^™ (76) 



fl2 T T I 



i?^" = (77) 

Substituting ([7^) -([771) into ( [75D we find the components of F^y. Further using the 
metric (|T3|) we find F^^ . We list here only those components of the (covariant and 
contravariant) curvature(s) that we will need below 





2CL 




c ■ 




+ 










r 


F- = 


2CL 




c ■ 


^2 


+ 


r . 



9 T 

(1 - CL)L' - — 

r 

9 T 

(1 + CL)U - — 

r . 



C^LL' 

X[iT.j]kXk — a;[jSj]o (78) 

C^LL' 

X[iX'j]kXk — £[iSj]o (79) 
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fcO 



-^—^kixi - -(1 - CL)L'S,o - C 



L" 



[l-CL)- 
r 



10 



^0) 



In addition to ( [78D ([79|) and (180|) , we will need the components F'*-''^, F^jj^ and Fmnro 
of the 4-form YM curvatures for the p = 2 case, which can readily be calculated using 
(^). These are needed for the calculation of the action density. 

Before giving the required action densities, we briefly verify that the fields given 
in ([781)- (pOD actally lead to self dual YM 2p-forms. The two selfduality equations can 
be stated HI as 



plijk 

Using the tensor-spinor identities 



kO 



__ lijkmnr rp 

" mnrO 



1 

3! 



^fcO 



^lijkmnr^mnrO 



^1) 

^3) 



in each case respectively, and (|78| 
equations and their solutions 



W), we find the following two simple differential 



L" 



2\ll 



2L 



I.e. 



1 

r 



I.e. 



3 

T 2 



^5) 
^6) 



for p = 1 and p = 2 respectively. These agree, through L = Iq, with (|T6D as expected. 

It follows from the selfduality equations (pi|)-(P^) that the action densities defined 
in (|38D and pO|) for the p = 1 and p = 2 cases simplify to 



5(1 



5, 



(1) 

(2) 



2 X Tr F'^Fij 
2 X Tr F^'^^F, 



i7) 



lijk 



in which the contributions of the terms 2Tt F^^Fko and 2Tt F"^"''^^ Fmnro are absorbed 

and (|88|). Thus for this purpose, we 



in the factors 2 on the right hand sides of (JS 

need only calculate the components F^\ Fij, F^^^^ and Fijki of the YM curvature, and 
not F^^^^ and FijkQ. In the definitions (|87|) and (^) we have omitted factors which 
would cancel with the angular volumes and the periods of the (Euclidean) time (^, 
which are necessary to render the C-P charges of the static spherically symmetric 
self-dual solutions equal to one. 
A direct calculation then yields 



5, 



(1) 



5, 



(2) 



2 X 2—{d 

I - 

2 X 2—{d-2){d 



2)^[2L'2 + 4(rf 



3) 



3)(rf-4)— -L^[L'^ + (rf-5) 



(89) 
(90) 
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where it is understood that d = A in ( |5UD and c? = 8 in (PPj), and we have left d in 
these expression by way of highhghting their relations to (^) and (^TJ). 

Substituting the solutions (^) and ( pUf ) in (pOf ) and (pO|), performing the angular 
integrations and the (Euclidean) time integrations over one period (P), there remains 
only to perform the radial integrations we to evaluate the action integrals 

C(p) J I 



S, 



(p) 



4) ^'""'dr 



(91) 



In (|9T|) , C(p) is a normalisation constant. The the angular integration over the d — 2 
dimensional angular volume gives 

27r('^-i)/2 

^'^'-'^ = m-l)/2) ^^^^ 

The factor P{p2), namely the period (^), is conributed by the (Euclidean) time integra- 
tion. C(p) is to be chosen such that the Chern-Pontryagin (C-P) charge of the selfdual 
spherically symmetric field configurations be normalised to unity. We calculate this 
normalisation factor for the Pi = P2 = 1 and pi = p2 = 2 cases in c? = 4 and d = 8 
dimensions respectively. The scale factor C in (^) and (0) is fixed to C = 1/2, as 
was done in the evaluation of the components of the Riemann tensor in (|73|)-(^). 

Substituting (^) for L{r) in (^), with d = 4, and performing the radial integral, 
we find 

C(i) = P(i)fi(2) . (93) 

Similarly performing the radial integral of ( pO]) with L{r) given by (^) and d = 8, 
we find 

C(2) = 90 P(2)^^(6) . (94) 

We now compare the values of the actions of the non-selfdual solutions by perform- 
ing the four and eight dimensional integrals of ( p9| ) and ( p] ) using the normalisations 
(p3|) and (@), respectively. The results of the actual integrations, for both one-node 
and two-node solutions respectively, are listed in ( |5B| ) and ( |5D| ) for the pi = 1 and 
p — 1 = 2 cases in that order. The resulting action integrals, analogous to (|9TD 



ip) 



with given by (|3; 



C(p) 

,iii), mM), 



5, 



(p) 



dr 



5 



(1) 

(2) 



0.959 
1.587 



(2) 
(1) 
(2) 
(2) " 



0(2 



take the values and (|93|),(|9 
= 0.992 
1.933 . 



(95) 

(96) 
(97) 



The superscripts in (p6|)-(p7|) pertain to the number of nodes, as in (pq),(p9|). 

The magnitudes of these actions are to be compared relative to the unit valued 
actions of their selfdual partners. In the latter case, these are simply equal to the unit 
C-P charges, while in the non-selfdual cases the C-P charges equal by virtue of the 
vanishing of FiQ, (^4]) . The only quantitative conclusion that can be made is, that it 
appears the actions grow with increasing number of nodes, which is not surprising. 
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IV Summary and discussion 



Since the known real solutions to the 4p-dimensional p-YM hierarchy are self-dual, 
it is of some interest to construct non-selfdual solutions, and this is what we have 
done by considering these systems on fixed curved backgrounds. YM fields on fixed 
curved backgrounds can be regarded as approximations to the fully interacting YM- 
gravitational fields, but in this paper we have restricted to the fixed gravitational 
backgrounds. Nonetheless we have discovered properties that are akin to those of 
the usual 4-dimensional Einstein-YM fields studied by Bartnik and McKinnon [ p^ . 
Specifically, we have found that there exist solutions for which the radial function 
parametrising the YM field has one- and two-nodes. It is quite likely that there 



should be a sequence of solutions, like in [jT3[, with increasing number of nodes, but 
we did not search for these. 

We have presented non-selfdual solutions to the pi-hierarchy of Yang-Mills (YM) 
systems on the fixed backgrounds of Schwarzchild and deSitter spaces. The specific 
constructions are made for the pi = 1 and the pi = 2 systems in various dimensions, 
but the qualitative conclusions arrived at are expected to remain valid for arbitrary 
Pi- 

The fixed Schwarzschild and deSitter metrics employed are the solutions of the 
vacuum gravitational equations of the p2-hierarchy of Einstein-Hilbert (EH) systems. 
Thus, for each of the two cases pi = 1, 2, we have employed p2 = 1, 2, 3 fixed back- 
grounds consistent with the requirement of asymptotic flatness of the metric. Since 
the deSitter metrics for all P2-EH systems are identical, this multiple choice of curved 
background is relevant only for the Schwarzschild metrics (^), given in But in 
that case (i.e. Schwarzschild) the non-selfdual solutions to the pi-YM systems we 
construct are real only in dimensions d = Api whence it follows that for the pi = 1 
case in 4-dimensions, only the p2 = 1 background can be used consistently with the 
asymptotic flatness condition. For the pi = 2 Schwarzschild case, we use all possible 
backgrounds p2 = 1,2,3. An interesting observation that can be made is, that for 
a given pi-YM system on a P2-EH background, the energy decreases with increasing 
P2 of the background. This was found to be true for both one-node and two-node 
solutions, cf. (|5|)-(|60|). 

In the pi = 1 Schwarzschild case, the one-node solutions are constructed both in 
closed form and numerically, while the two-node solutions of this system, as well as 
all (one- and two-node) solutions in the p = 2 case were possible to construct only 
numerically. The profiles of the functions K{r) are exhibited in Figs. 1,2. 

In the pi = 1 deSitter case, the one-node solutions are constructed both in closed 
form and numerically, while the solutions in the pi = 2 deSitter case these are con- 
structed only numerically. No solutions with more than one mode were found in these 
cases numerically, but we do not know if this indicates the non-existence of such so- 
lutions. If they exist, then it would be a challenge to find these, but it is possible that 
they do not since the two backgrounds, Schwarzschild on which two-node solutions 
exist and deSitter, are qualitatively quite different from each other. Notably the in- 
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tervals on which they are defined are respectively non-compact and compact. In the 
deSitter cases we did not exhibit the profiles of the functions K{r) for these solutions 
as these all have one-node only, and, they are very close to the closed form solution 
of the pi = 1 case. Instead we have listed their properties, namely the position of the 
(single) node and the energy integral, in Tables 1 and 2, for the pi = 1 and pi = 2 
cases respectively. 

That the solutions discussed above are non-selfdual is not a matter of note un- 
til one considers that in the specific dimensions 4pi = 4p the pi = p-YM fields on 
double-selfdual (2p-form Riemann tensor) p2 = p-EH background, are (single) selfd- 
ual. In these cases where pi = P2 = P, the action of the gravitational system vanishes 
due to the vanishing of the stress tensor of the selfdual YM fields, so that the to- 
tal action equals the Chern-Pontryagin (C-P) charge. In the spherically symmetric 
Schwarzschild cases considered, this is the unit C-P charge. The comparison of this 
action with the action pertaining to the corresponding non-selfdual solution reveals 
an interesting feature. (Note that the non-selfdual solutions have zero C-P charge.) 
It was found that for the 4-dimensional [pi = 1) YM system on p2 = 1 Schwarzschild 
background, the actions of the nonselfdual solutions with one and two nodes respec- 
tively were equal to 0.959 and 0992, which are slightly less than the unit selfdual 
action. (Of these the exact value of the one-node action was found in Ref. [|1|.) We 
have found that for the pi = 2 YM system on p2 = 2 Schwarzschild background in 
8 dimensions, the actions of the one node and two node solutions respectively are 
equal to 1.587 and 1.933, which are appreciably larger than the unit magnitudes of 
the corresponding selfdual actions. 

In addition, we have exhibited anti-deSitter solutions of the p-YM systems in 4p di- 
mensions which can be related with meron-type solutions in fiat space-time through 
conformal transformations. Very recently, soliton-types of solutions have been ob- 
tained for the Einstein- Yang-Mills equations with asymptotically anti-de Sitter space 


There remain two obvious directions in which the study of the present work can 
be developed to complete and enhance the conclusions drawn here. The most obvious 
extension is to find the (non-selfdual) solutions of the full YM-EH systems taking into 
account fully the backreaction of gravity on the gauge fields, rather than studying the 
gauge fields on a fixed background. The actions of the solutions in that case would 
differ from what we have found, but maybe not very much. What is important is that 
the topological (C-P) charges of these solutions will still be equal to zero so that like 
those on the fixed backgrounds, these solutions will also describe (unstable) saddle 
points. Another direction, is to study the solutions of combined p2 EH systems in 
dimensions d > 4p2. (Note that in d = 4p2 the P2-EH system (|^)-(|^) reduces to the 
Euler-Hirzebruch topological density which has no dynamics.) It would be interesting 
to find out how the inclusion of a higher p2-EH term (presumably with a small coupling 
constant) would modify the metric and other properties of the leading p2 system. In 
addition, one can contemplate the addition to these types of systems, members of the 
Pi-YM hierarchy. 
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Appendix: A nonselfdual solution for AdS^ back- 
ground. 

For (i = 4 a divergent solution in deSitter background was shown |l|, |^ to be related 
through conformal transformations to meron-type solutions in flat space-time. This, 
particularly simple, solution corresponds (with A > 0) to 

N={l-Ar'^), K = N-^ = {I - Ar^)-^ (A.l) 

Here we note that changing the sign before A one obtains the anti-deSitter case {AdS^) 
and 

N={l + Ar^), = = (1 + Ar2)-5 . (A.2) 

This provides a solution of ( ^6l) with p = 1 and d = 4, namely of 

{NKJ = r-^K{K^ - 1) (A.3) 

where now K is no longer singular. 

Remarkably, we note that ( |A.2[ ) is a solution to the p-YM system in 4p dimensions, 
satisfying the equation 

[N{K'' - ly-^K']' = {2p - l)r-^K{K^ - Vf (A.4) 

that results from (|36|) setting c/ = 4p in it. This is not surprising since it is known 
that the p-YM system on R^^ supports meron solutions . 

One can now evaluate the total actions for these 4p-dimensional solutions (|A.2| ) 
and study their different properties. We will restrict this to the (j9 = 1) 4-dimensional 
case, setting A = 1 henceforth in this Appendix for simplicity. 

The radial integral of (^), in the absence of the horizon, should now be replaced 

by 

/(4) = fdxx-'{2N{x^y + {k' - if) (A.5) 



where 

One obtains 



N = {l + x^), K={l + x^)-^ 
^4) = 3 / ^,,, dx = ^ (A.6) 



[l + x^)^'" 4 
Thus, corresponding to (^9]), one obtains a finite spatial integral 

37r2 



(A.7) 

The factor from the time integration depends on the chosen context. Now there is no 
horizon to be desingularized and the discussion of Sec.l is not directly relevant. But 
one can start by considering the hypersurface 

-tl-tl + xl + xl + xl = -l (A.8) 
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In terms of the spherical coordinates 



{hM) ^{T,^) (A.9) 

the metric on the hypersurface 

^2_rp2^ _1 

is 

ds'^ = -(1 + r2)#2 + (1 + r2)-Mr2 + r^rfl^s (A.IO) 
In this context the ^/^-integration gives a factor 2tt and one obtains a total action 

3tt^ (A. 11) 

But often it is preferable to consider the covering space (CAdS) replacing ip E hy 
t e R. Then the action is evidently divergent. 
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Figure Captions 

Figure 1 The profiles of two solutions of Eq. (|37|) (one with one node and one with two 
nodes) for the case pi = l,p2 = I. 

Figure 2 The profiles of two solutions of Eq. (|37|) (one with one node and one with two 
nodes) for pi = 2 and for P2 = 1 (solid), P2 = 2 (dashed), P2 = 3 (dotted). 
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